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Abstract. Given a fc-uniform hypergraph on n vertices, partitioned in k equal parts 
such that every hyperedge includes one vertex from each part, the fc-Dimensional Match- 
ing problem asks whether there is a disjoint collection of the hyperedges which covers 
all vertices. We show it can be solved by a randomized polynomial space algorithm in 
0*(2 n(fc_2)/fc ) time. The 0*() notation hides factors polynomial in n and k. 

The general Exact Cover by fc-Sets problem asks the same when the partition constraint 
is dropped and arbitrary hyperedges of cardinality k are permitted. We show it can be 
solved by a randomized polynomial space algorithm in 0*(Cfc) time, where C3 = 1.496, C4 = 
1.642, C5 = 1.721, and provide a general bound for larger k. 

Both results substantially improve on the previous best algorithms for these problems, 
especially for small k. They follow from the new observation that Lovasz' perfect matching 
detection via determinants (Lovasz, 1979) admits an embedding in the recently proposed 
inclusion-exclusion counting scheme for set covers, despite its inability to count the perfect 
matchings. 



1. Introduction 

The Exact Cover by fc-Sets problem (XfcC) and its constrained variant fc-Dimensional 
Matching (A;DM) are two well-known NP-hard problems. They ask, given a /c-uniform hy- 
pergraph, if there is a subset of the hyperedges which cover the vertices without overlapping 
each other. In the /cDM problem the vertices are further partitioned in k equal parts and 
the hyperedges each includes exactly one vertex from each part. While being two of the 
21 items of Karp's classic list of NP-complete problems [6] for k > 3, little is known on 
their algorithmic side. In this paper, we present stronger worst case time bounds for these 
problems by combining Lovasz' perfect matching detection algorithm via determinants [TO] 
with the inclusion-exclusion counting for set covers pp. We show 

Theorem 1.1. k-Dimensional Matching on n vertices can be solved by a Monte Carlo 
algorithm with exponentially low probability of failure in n, using space polynomial in n, 
running in 0*(2 n ^ k ~ 2 ^ k ) time. 

Theorem 1.2. Exact Cover by k-Sets on n vertices can be solved by a Monte Carlo al- 
gorithm with exponentially low probability of failure in n, using space polynomial in n, 
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Table 1: Comparison of the base c in the 0*(c n ) runtime of previous and the new algo- 
rithms. 

running in 0*(c%) time, with C3 = 1.496, C4 = 1.642, C5 = 1.721, c§ = 1.771, c-j = 1.806, and 
in general c k < 2 (8A15k°- g - k (k - l) os (k - l^)* -1 - 5 ) -17 * 

These bounds are large improvements over the previously known ones. In particular, 
for three dimensional matching our algorithm runs in time asymptotically proportional to 
the square root of the previous best algorithm's runtime. 

We hope the present paper conveys the message that inclusion-exclusion is amendable 
not only to counting problems, but can at times be used more directly to settle the decision 
version of a problem. 

1.1. Previous Work 

Perhaps the most famous algorithmic contribution on the subject of exact covers is 
Knuth's Dancing Links paper |7J, which actually just addresses a general implementation 
issue which saves a small constant factor in the natural backtracking algorithm for the 
problem. About the backtracking approach on exact cover he writes "Indeed, I can't think of 
any other reasonable way to do the job in general" . While we certainly may agree depending 
on how much you put in the words "reasonable" and "general" , we must point out that the 
best provable worst case bounds for the problems are obtained by analyzing very different 
algorithms. Bjorklund et al. [2] uses inclusion-exclusion and fast zeta transforms on the full 
subset lattice to show that exact set covers of any n vertex hypergraph can be counted in 
O* (2 n ) time even when the number of hyperedges to choose from are exponential. Restricted 
to /c-uniform hypergraphs, Koivisto [8] proposes a simple clever dynamic programming over 
subsets which show that Exact Cover by fc-Sets can be solved in o*(2 n ( 2fc - 2 )/v / (2fc-i) 2 -2ln(2)) 
time. The algorithm is actually capable of counting the solutions and also works for not 
necessarily disjoint covers. It does, however, use exponential space. The best previous 
algorithm for the problem using only polynomial space is given in [1] and has a runtime 
bound in 0*((l+k/(k— l))™( fc ~ 1 )/ fc ). For /c-Dimensional Matching, the best known algorithm 
as far as we know is an 0*(2 n ( k ~ 1 ^ k ) time algorithm resulting from a generalization of 
Ryser's inclusion-exclusion counting formula for the permanent |12| . presented in [1J. A 
comparison of the bounds guaranteed by these algorithms and the ones given in this paper 
is shown in Table [1] for small k. 

For k = 2 the problems X2C and 2DM are better known as the problems of finding a 
perfect matching in a general and bipartite graph, respectively. For these problems several 
polynomial time algorithms are known. We definitely admit that it seems like an obvious 
idea to try to reduce the k > 2 cases to the k = 2 case searching for faster algorithms 
for larger k. Still, we believe that it is far from clear how to achieve this efficiently. In 



EXACT COVERS VIA DETERMINANTS 



97 



this paper we make such an attempt by reducing the k > 2 cases to variants of one of the 
first polynomial time algorithms for detecting the existence of perfect matchings: Lovasz' 
algorithm from [TO] which evaluates the determinant of the graph's Tutte matrix [13J at a 
random point. 

2. Our Approach 

2.1. Preliminaries 

We use the terminology of (multi)hypergraphs. A hypergraph H = (V, E) is a set V of 
n vertices and a multiset E of (hyper) edges which are subsets of V . Note in particular that 
with this definition edges may include only one (or even no) vertex and may appear more 
than once. In a fc-uniform hypergraph each edge e S E has size |e| = k. Given a vertex 
subset U C V, the projected hypergraph of H = (V, E) on U, denoted H[U] = (U, E[U]) is a 
hypergraph on U where there is one edge ejj in E[U] for every e S E, defined by ejj = efl U, 
i.e. the projection of e on U. 

We study two related problems. 

Definition 2.1 (fc-Dimensional Matching, fcDM). 

Input: A /c-uniform hypergraph H = (Vi U V2 U • • • Vk, E), with E C Vi X V2 X • • • Vk- 
Question: Is there SCB s.t. U se5 s = U V% U • • • Vk and Vsi / s 2 G 5 : si n s 2 = 0. 

Definition 2.2 (Exact Cover by fc-Sets, XkC). 
Input: A /c-uniform hypergraph H = (V,E). 

Question: Is there S C E s.t. L^ggs = V and Vsi 7^ s 2 G 5* : si H s 2 = 0. 
For a matrix A we will by Ajj denote the entry at row i and column j. 

2.2. Determinants 

The determinant of an n x n-matrix A over an arbitrary ring R can be defined by the 
Leibniz formula 

n 

det(A)= sgn(a)l[A iMi) (2.1) 

er:[n]— >[n] 1=1 

where the summation is over all permutations of n elements, and sgn is a function called 
the sign of the permutation which assigns either one or minus one to a permutation. In this 
paper we will restrict ourselves to computing determinants over fields of characteristic two, 
GF(2 m ) for some positive integer m. In such fields every element serves as its own additive 
inverse, and in particular so does the element one, and the sgn function identically maps 
one to every permutation. Thus it vanishes from Eq. 12.11 in this case, and the determinant 
coincides with another matrix quantity, called the permanent: 

n 

per(A)= II A v^) (2.2) 

<r:[n]-»|n] i=1 

Permanents of 0-1-matrices over the natural numbers are known to count the perfect match- 
ings of the bipartite graph described by the matrix. The reader may subsequently be 
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tempted to think that this identity of determinants and permanents over fields of char- 
acteristic two is the property that makes our algorithms work. There is however nothing 
magical about these fields in this context. Our reason for working in GF(2 m ) is simply that 
with this choice of fields we don't even have to define the sign function, making several of 
the proof arguments later on much easier to digest. In principle though, any large enough 
field will work, with slightly more complicated proofs. 

The interesting property of the determinant that we will exploit here is that although 
it is defined above in Eq. 12.11 as a sum of an exponential number of terms, it admits 
computation in time polynomial in n. This can be achieved for instance via the so called 
LU-factorization of the matrix which almost any textbook on linear algebra will tell you. 
In fact, computing the determinant is no harder than square matrix multiplication, see [3], 
and hence it can be done in 0(n w ) field operations where oj = 2.376 is the Coppersmith- 
Winograd exponent 

2.3. Inclusion Exclusion for Set Covers 

Let us review the inclusion-exclusion counting scheme for exact set covers presented by 
Bjorklund and Husfeldt in [lj: Given a /c-uniform hypergraph H = (V,E) and any subset 
U C V, we can count the number of Exact Covers by fc-Sets, denoted #X.kC(H), by the 
inclusion-exclusion formula 

#XkC(H)= Yl (-l) lXl W(H,U,X) (2.3) 

XQV-U 

where W{H, U, X) counts the number of ways to exactly cover U with | V|//c edges in H[U] 
whose corresponding edges in H are disjoint from X. Put differently, W(H,U,X) counts 
the number of ways to pick |V|/A; edges from H, all having an empty intersection with 
X, which cover U without any overlap. In particular, when U = it is straightforward to 
compute W(H, 0, X) by just counting the number of edges in H disjoint from X, calling this 
quantity d(X), and then computing the binomial ( d ^)- In pQ, some examples where this 
algorithm could be accelerated by choosing a larger U were identified where the speedup 
was obtained by utilizing U's such that the projected hypergraph on U had low path-width. 
This enabled efficient counting by dynamic programming over a path decomposition. 

2.4. Moving to GF(2 m ) 

In this paper, we find a new way to allow a large U to expedite the computation of the 
formula Eq. 12.31 above. We observe that whenever the projected hypergraph contains edges 
of size at most two, we can use determinants to compute the formula faster. We note that 
if the problem of counting perfect matching had an efficient algorithm A, we would almost 
immediately get an 0*(2 n ( k ~ 2 ^ k ) time algorithm for the /cDM problem. We would simply 
let U be any two of the parts in the input partition, and use A to compute W(H, U,X). 
Unfortunately, counting perfect matchings even in bipartite graphs is ^P-complete |14j . 

The key insight of the present paper circumvents the apparent obstacle formed by the 
intractability of counting matchings: we only need to be able to efficiently compute some 
fixed weighted sum of the matchings (with no weights set to zero) . This is exactly where the 
determinants come to our rescue. The price we pay is that we have to give up counting the 
solutions over the natural numbers. Here we demonstrate the result through counting over 
fields of characteristic two which only allow us to detect if there is a cover at all and gives us 
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little knowledge of their number. Furthermore, to avoid having an even number of solutions 
cancel we will employ a fingerprint technique, very much in the same spirit as Williams [15] 
recently extended the fc-path detection algorithm based on an algebraic sieving method of 
Koutis [9]. The fingerprint idea is to think of the computation as evaluating a polynomial 
of a degree much smaller than the number of elements of its base field and then computing 
it at a randomly chosen point. The fact that a polynomial cannot have more roots than 
its degree assure us that with great probability we discover with this single point probing 
whether the polynomial is the zero-polynomial or not. We will in fact use the multivariate 
polynomial analogue, see e.g. [TT] . 

Lemma 2.3 (Schwartz-Zippel). Let P(xi,x 2 , x n ) be a non-zero n-variate polynomial of 
degree d over a field F. Pick r\,r 2 , ■~,T n £ F uniformly at random, then 

Pr(P(n,r 2 ,...,r n ) = 0) < — 

For now, it is sufficient to think of the inclusion-exclusion formula of Eq. 12.31 as evalu- 
ating a multivariate polynomial over the base field GF(2 m ) for some m. In what follows we 
will associate with all edges e in the input hypergraph a variable v e . Our modified version 
of Eq. 12.31 reads as follows. 

Lemma 2.4. Given an XkC-instance H = (V,E) and the family of all its solutions S, we 
have that, for every subset U C V , 

£ W 2J (H,U,X)= £ H^ (e) (2-4) 

XCV-U E'^Se&E' 

where the computation is over a multivariate polynomial ring over GF(2 m ), f is a function 
mapping the edges to the positive integers, and 

W 2J (H,U,X) = Y^ II v e {e) ( 2 - 5 ) 

E" ee£" 

where the summation is over all E" C E, satisfying four constraints 

• Avoidance, Ve G E" : e n X = 

• Cardinality, \E"\ = \V\/k 

• Coverage, U C U eg £//e 

• Disjointness, Vei ^ e 2 6 E" : e\ n e 2 (1 U = 

Proof. First, note that every E' G S fulfills all four conditions Avoidance, Cardinality, 
Coverage, and Disjointness for X = 0, but violates Avoidance for every other X, irrespective 
of the choice of U. Thus, the contribution HeeB' Ve °^ ever Y solution E' is counted precisely 
once. 

Second, a non-solution E" obeying the three conditions Cardinality, Coverage, and 
Disjointness, fulfills the Avoidance condition for an even number of choices of X irrespective 
of U, namely for all subsets of the elements of V that the union of the sets in E" fails to 
cover. Hence, all of these contributions n e e£" v ^ cance l each other since we are working 
in a field of characteristic two. ■ 



100 



A. BJORKLUND 



Combining the two Lemmas above 12,31 and 12,41 into an algorithm choosing a random 
point ri,T2, ■•■j^ibi £ GF(2 m ) and evaluating the left-hand sum of Eq. 12.41 in the straight- 
forward fashion, we get: 

Corollary 2.5. Given an XkC-instance H = (V, E) and a subset U C V , there is a Monte 
Carlo algorithm which returns "No" whenever there is no cover and returns "Yes" with 
probability at least 1 — max eg E f(e)\V\/(k2 rn ) when there exists at least one, running in 
time 0*{2\ v \~\ u W(W2j, U)), where r(W2j,U) is the time required to evaluate any of the 
polynomials W2j(H, U,X) for X C V — U, in a random point over the base field GF(2 m ). 

Note that by letting m be in the order of n, when / is bounded by a constant, we 
get exponentially low probability of failure in n. Armed with Corollary 12.51 we can start 
looking for projections U over which the computation of W2j(H,U, X) is easy. The next 
two sections will describe two examples of how we can use determinants to accelerate the 
computation. 

3. /c-Dimensional Matching 

We begin by the easier application, A:DM. For this problem we can trivially find a large 
vertex subset on which the projected instance is a multigraph, and in fact also bipartite: we 
just use any two of the parts in the vertex partition given as input. Edmonds [5j observed 
that one could relate a bipartite graphs' perfect matchings to the determinant of a symbolic 
matrix. A perfect matching is a collection of disjoint edges so that every vertex is covered by 
precisely one edge. To a given a bipartite graph G = (U, V, E), n = \U\ = \V\, he associated 
an n x n-matrix A with rows representing vertices in U, and columns the vertices of V, 
and equated Ajj with a variable Vij if (i,j) £ E and zero otherwise. He showed that the 
determinant of A is non-zero iff G has a perfect matching. We will use essentially the same 
result, with the small exception that we need to deal with multiple edges between a vertex 
pair, making sure all contributes. Formally 

Definition 3.1. Given a hypergraph H = (V,E) and a subset U C. V such that the 
projected hypergraph H[U] is a bipartite multigraph on two equally sized vertex parts 
U\ U U% = U, its Edmonds matrix, denoted E(ff, U%, U2), is defined by 

E(H,U 1 ,U 2 ) l , 3 = J2 v - 

e=(i,j)eE[U] 

ieu u jeu 2 

where again, v e is a variable associated with the edge e. 

We formulate our Lemma in terms of a special case of XfeC instead of kDM directly to 
capture a more general case. 

Lemma 3.2. For a XkC-instance H = (V,E) and two equally sized disjoint vertex subsets 
U\, U2 Q V such that the projected hypergraph H[U\ U U2] is a bipartite multigraph, 

det(B(H,U 1 ,U 2 )) = ]T 11 v e (3.1) 

where the computation is over a multivariate polynomial ring over GF(2 m ) for some m and 
the summation is over all perfect matchings M. in H\U\ U U^]. 
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Proof. By definition of the determinant 12,11 the summation is over all products of n of the 
matrix elements in which every row and column are used exactly once. Transferred to the 
associated bipartite graph, this corresponds to a perfect matching in the graph since rows 
and columns represent the two vertex sets respectively. Moreover, the converse is also true, 
i.e. for every perfect matching there is a permutation describing it. Hence the mapping is 
one-to-one. The inner product counts all choices of edges producing a matching described 
by a permutation a since: 

n n 

ri E ^)=n E «"= e U V * ( 3 - 2 ) 

i=l i=l e=(i,<r(t)) MeM(cr) e£M 

where JVi(cr) is the set of all perfect matchings e±, 62, e n such that a = (i, o~(i)). m 

3.1. The Algorithm 

Now we are ready to prove Theorem 1 1,11 Given an input instance H = (Vi, V%, V^, E) 
to the A:DM problem where Vi, V2, describe the vertex partition of the n vertices, we 
simply let U = V\ U V2 in the algorithm described by Corollary 12.51 with / mapping one 
to every edge. To compute W2j(H, U,X) we construct the Edmonds matrix of the hyper- 
graph H restricted to its edges disjoint to X, projected on U, with the variables replaced 
by the random sample point (n, r%, t\e\) chosen. Next we compute its determinant. 
The correctness follows from Lemma 13.21 after noting that every perfect matching in a 
projected hypergraph contains n/k disjoint edges. The runtime bound is easily seen to be 
*( 2 n(fc-2)/fc) gince = 1^ + |y 2 | = 2n /k. 

4. Exact Cover by A;-Sets 

Next we proceed to the X£;C problem. In comparison to the fcDM we are faced with a 
number of additional obstacles on our way to a similar result. 

• First, a projection will typically capture edges differently, some will have large pro- 
jections and some no at all. 

• Second, in particular the projected edges will probably not form a multigraph. 

• Third, even if they did it may not be a bipartite one. 

For the first obstacle, we will prove that it is sufficient to find a projection on which 
at least one cover's edges all leave projected edges of size two or less. This is basically 
an extension of the idea for the X/cC algorithm in proposition 10 in |lj. There, a vertex 
subset U is picked uniformly at random of a carefully chosen size, and in the projected 
hypergraph only the edges which leave a projection of size one or less are kept. Then the 
inclusion-exclusion formula Eq. 12.31 is used after noting that W(H,U,X) is now easy to 
compute. The process is repeated a number of times dictated by the size of U. The best 
size to use is a trade-off of the resulting summation runtime and the probability that a cover 
is projected gracefully in the sense that all its edges are kept after the projection. 

For the second obstacle, in addition to handling multiple edges we also need to count 
perfect matchings in which loops, i.e. edges connecting a vertex to itself, count as covering 
the vertex of its endpoints. Since this means that not all perfect matchings will involve 
the same number of edges, we have to take special care to make the determinants useful. 
We use polynomial interpolation to solve for the contributions of matchings of the same 
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size separately to be able to fulfill the Cardinality constraint for Wij in Corollary 12.51 To 
this end we introduce an auxiliary variable s parametrizing the matrices and use several 
determinant calculations. 

For the third obstacle, we will use a variation of a result generalizing Edmonds' due to 
Tutte [13]. He showed that even for general not necessarily bipartite graphs one can make 
a connection between its perfect matchings and the determinant of a symbolic matrix, 
although twice as large matrices in both directions are required. To a given a graph G = 
(V,E),n = \V\ he associates an n x n-matrix A with rows and columns representing the 
vertices, and assigns Aj = Vij for i < j and Ajj = — Vij for i > j with Vij a variable 
for each edge G E. The remaining entries are set to zero. The determinant of A is 
non-zero iff G has a perfect matching. 

We define matrices similar to Tutte's: 

Definition 4.1. Given a hypergraph H = (V,E) and a subset U C V such that in the 
projected hypergraph H[U] all edges have size at most two, its Tutte matrix of index s, 
denoted T^ S \H, U), is defined by 



Lemma 4.2. For a XkC-instance H = (V,E) and a vertex subset U C V such that in the 
projected hypergraph H[U], every edge has size at most two, 



where the computation is over a multivariate polynomial ring over GF(2 m ) for some m, the 
summation is over all perfect matchings Ai in H[U], A(M) is the number of loops in the 
matching M , and p{e) = 1 if e is a loop and p{e) = 2 otherwise. 

Proof. By definition of the determinant 1 2 . 1\ the summation is over all products Y\7=i "^i 
for a permutation a. Call a permutation a good if Vi : a{a{i)) = i holds, and bad otherwise. 
We will argue that only good permutations contribute to the sum. To see why, consider a 
bad a. Then there exists a smallest i such that a(a(i)) = j 7^ i. Look at the cyclic sequence 
{q} where cq = i and c& + i = er(cfc) for k > 0. Let L > 2 be the smallest positive integer 
such that cl = i (Note that there must be one and that all q in between must by distinct 
since every element in 1 through n is mapped to exactly once). Next define a cycle reversal 
operation D mapping bad permutations on bad permutations by letting D(a) be identical to 
a except in the points c\ through cl, where instead D(a)(ci) = Ci-\. Now first observe that 
the reversal operation is dual in the sense that D(D(a)) = a and that D(a) 7^ a since L > 2, 
and hence every bad permutation can be uniquely paired with another bad permutation. 
Second note that the contribution of a bad permutation is identical to the contribution of 
its dual, since the Tutte matrices are symmetrical. Thus, since we are counting in a field of 
characteristic two, they cancel each other. 

Next we continue to observe that the good permutations describe precisely the structure 
of all possible perfect matchings in a multigraph: z's such that a{i) 7^ i describe ordinary 
two-vertex edges in the matching, and i's such that o~(i) = i describe loops. 





(4.1) 
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The inner product of Eq. 12.11 reads 

n T S<.,=( n •!>)( n e «.)= e » A(M) n»« <«> 

i=l \i,i=<r(i) e=(i,i) / \i,i^<j(i) e=(i,<r(i)) / AfeM(ir) e£M 

where .M(o") is the set of all directed perfect matchings e%,e2, ■ ■■,e n described by the good 
permutation a for which ei = (i,o~(i)). 

Now consider a directed perfect matching e%, ei, e n such that for some j, ej ^ e <r(i)i 
and refer to it as being bad. We will see that all of these cancel in very much the same way 
as the bad permutations did. Namely, again find the smallest j for which this is the case, 
and define a reversal operation R a mapping bad directed perfect matchings onto themselves 
by exchanging ej and e-aU) ■ Since this operation pairs up the bad directed perfect matchings 
(i? CT ({ej}) ^ {e^ and R a (R a ({ei})) = {ej}) and we work in a field of characteristic two, 
their contributions cancel. Thus we are left with only good permutations and good directed 
perfect matchings. The latter can be thought of as undirected perfect matchings in which 
every non-loop edge is included twice in the product. ■ 

To find the contributions of matchings of the same size separately, think of the match- 
ings partitioned in groups Mo, Mi, ...,M n according to the number of loops of the match- 
ing. We can rewrite the determinant in Lemma 14.21 as 

n 

det(T« (H, U)) = J2 siM i ( 4 - 3 ) 
i=0 

where Mj = YImeM- TieeM u e are ^ ne quantities we seek. The right hand side of Eq. 14.31 
is a degree n polynomial in s and thus we can solve for Mq,M\, ...,M n by computing 
det(TW(fT, U)) in n different choices of s, and use Lagrange's interpolation formula to 
recover the sought values. In fact, either there are no matchings with an odd number 
of loops or no matchings with an even number of loops depending on the parity of \U\. 
Consequently, the evaluation of n/2 points suffices, but we disregard from this optimization 
possibility for simplicity. Once we have the Mi's we are close to be able to compute W2J 
efficiently according to the following Lemma: 

Lemma 4.3. Given a XkC instance H = (V,E) and a U C V such that for all edges 
e £ E,\e Pi U\ < 2, /(e) = 2 if \e f] U\ =2 and 1 otherwise, and any X C V — U , 

W 2J (H,U,X) = i>(^ - l^^m (4.4) 

where Mi = ^MeA4 TleeAf v ^ are ^ e contribution of all matchings Mi containing 
exactly i loops in the projected hypergraph of H on U restricted to the edges disjoint to X, 
and 

z ® = e n ^ ( 4 - 5 ) 

E'CZ e'£E> 
\E'\=i 

with Z defining the set of edges e disjoint to X also having an empty intersection with U. 

Proof. The Mj's count the contribution of all ways to cover U with the edges which leaves 
a non-empty projection on U and the Z(z)'s count the contribution of all ways to choose 
edges leaving an empty projection. Note that a matching from Mi involves exactly [ ^ j 
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edges if it exists. The right hand side of Eq. 14.41 convolutes over all ways their total number 
of edges could equal |V|/fc in order to meet the Cardinality constraint in Corollary 12.51 ■ 

The only piece missing is a simple way to evaluate Z(i), and we note that it can be 
done by dynamic programming through a simple recursion. Number the edges in Z defined 
in Lemma [431 arbitrarily as ei, e2, e p , set Zi = {ei, e%, ej}, and define 



z(i,j) 



E 

E'c: 
\E'\- 



n 



These can be solved for by 

1 


z(i - 1, j - l)v ej + z(i,j - 1) 
and we finally compute Z(i) through Z(i) = z(i,p). 



i = j = 
i = or j 
otherwise 



(4.6) 



(4.7) 



4.1. The Algorithm 

We are ready to prove Theorem 11.21 First we describe the algorithm. Given an input 
instance H = (V, E) to the X/cC problem, we compute two parameters t and I depending 
on k. These are given by the calculations in the next section [4~2l We repeat the following 
procedure until we detect a cover, in which case we report so, or have tried unsuccessfully 
i" times, in which case we report that no cover was found: 

Algorithm 4.4. 

(1) Choose a in-sized subset U QV uniformly at random. 

(2) Construct H v = {V, E v ) where E u = {e\e G E, \e n U\ < 2}. 

(3) Run the summation algorithm in Corollary 12.51 on Hjj, using U, and let /(e) = 2 if 
|e n U\ =2 and 1 otherwise. Use the method of the previous section [3] to compute 
W 2J (Hu,U,X), i.e. 

(a) Construct the Tutte matrices of Hjj[U] restricted to its edges which are 
disjoint to X for s = g l , < i < \U\ where g is a generator of the multiplicative 
group in GF(2 m ). 

(b) Compute the determinants of T( s ). 

(c) Use Lagrange interpolation to solve for the Mj's via Eq. 14.31 

(d) Calculate the Z(i)'a by Ea.K7\ 

(e) Evaluate Eq. Ol 

Given that the random U is such that all edges in some exact cover S are kept in 
Hjj, the previous Section U] verifies its correctness: Lemmas 14.21 and 14.31 together with the 
observation that g l for 1 < i < \U\ are all distinct points, assures us that step (3) of the 
algorithm works. We are left with deciding t and / to make it very likely that some exact 
solution is kept at least once and tune them to get the best possible runtime. 



EXACT COVERS VIA DETERMINANTS 
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k 


712 


72 


t 


jl/n 


Cfc 


3 


0.961 


0.679 


0.547 


1.092 


1.496 


4 


0.936 


0.613 


0.387 


1.073 


1.642 


5 


0.921 


0.583 


0.301 


1.060 


1.721 


6 


0.912 


0.565 


0.246 


1.050 


1.771 


7 


0.905 


0.554 


0.208 


1.043 


1.806 


8 


0.900 


0.546 


0.181 


1.038 


1.832 



Table 2: Numerically found parameters t\ 2 and T2 which approximately minimizes Cfc. 



4.2. Runtime Analysis 

Our runtime analysis hinges on the probability that any fixed solution S to the XfcC 
instance H = (V,E) when projected on a subset U C V chosen uniformly at random from 
the tn-sized subsets of V for some fraction t of the vertices, gets all its edges to leave a 
small projection on U, namely Ve S S : |e D U\ < 2. We denote this event by e(t). If 
we repeat the process I times, the probability that none of the I independent random 
selections for U is successful in the sense that they retain S after the projection, is at 
most (1 — Pr(e(t))Y < e~ ^■ £ ^ 1 . Consequently, we need / = log(e _1 )Pr(e(t)) _1 to get 
probability at least 1 — e for one or more of the I selections to be successful. Thus we may 
use e = c~ n for some constant c > 1 to get an exponentially low probability in n of failure 
without increasing the number of repetitions / by more than a polynomial factor. 

To bound the probability of the event, we count the number of good tn-sized subsets 
of the vertices. This is a binomial sum (actually a trinomial one) over the number of edges 
in the solution S which gets a projection of size two: 

To lower bound this sum of all non-negative terms, we will use just one of them. Let 
N = n/k and parametrize tn = T\ 2 N + t 2 N where T12 is the fraction of sets in the solution 
S which gets at least one of its elements chosen, and T2 is the fraction of sets that gets two. 

Then, we bound our probability as the quotient of the single term lower bound on the 
number of good sets and the number of all sets ( t ™) to 



2r 2 N ( kN 



\(t 12 +t2)n) 

The runtime of Corollary 12 .51 is 0*(2 n ~ tn ) given our polynomial time algorithm for com- 
puting W2j- Omitting polynomial factors, Algorithm 14.41 for XA;C has to run for Pr(e(t)) _1 
different choices of U in the worst case. Let Tf. t denote the final runtime, and expand the 
binomials of Eq. 14.91 to get: 

2 n-tn 2 kN -^ N (T 2 N)\(N - T 12 N)\(T 12 N - T 2 N)\(kN)\ 

M - Pr(e(t)) ~ N\k T ^ N {k - l) T2N (kN - t 12 N - t 2 N)\(t 12 N + t 2 N)\ [ ' ' 



If we replace the factorials with Stirling's approximation n! G 8(y/n(n/e) n ) and divide 
(N/e) k+1 out of both numerator and denominator, we are left with a slightly less intimi- 
dating expression 
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T < ( 2( fc -^)r 2 r2 (r 12 - r 2 )^-^(l - ria) 1 -^ 

k,t ~ I k T ^~ k (k - 1) T2 (A; - T12 - T 2 ) fe " ri2 " T2 (n2 + T2) ri2+T2 

Rewriting this as T^^ < cJJ we see that can be obtained as the k:ih root of the 
expression within the brackets in Eq. 14.111 Solving numerically for the choices of t\% and t 2 
that minimizes Ck we find that the minimum moves slightly with increasing k, see Table EJ 
The minimum, however, lies in a quite flat neighborhood within a large vicinity of the actual 
minimum, and comparable bounds not too far from the best possible with our technique are 
obtained with fixed parameters for all k by, say, t\% = 0.9 and r 2 = 0.6. With this choice of 
parameters in Eq. 14.111 we obtain the general bound in Theorem 11.21 
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